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The present invention relates to methods for treatment of digitised audio 
signals (digital stored sample values from an analogue audio waveform signal) 
and, in particular (although not exclusively) to the application of such methods to 
extending the duration of signals during playback whilst maintaining or modifying 
their original pitch. The present invention further relates to digital signal 
processing apparatus employing such methods. 

The enormous increase in multimedia technologies and consumer 
expectation for continually higher standards from home audio and video 
systems has led to a growth in the number of features available on home 
multimedia products. These features are vital for product differentiation in an 
area that is extremely cost sensitive, and so new features are usually 
constrained with critical CPU and memory requirements. 

One such feature is slow motion audio based around a Time Scale 
Modification (TSM) algorithm that stretches the time content of an audio signal 
without altering its spectral (or pitch) content. Time scaling algorithms can either 
increase or decrease the duration of the signal for a given playback rate. They 
have application in areas such as digital video, where slow motion video can be 
enhanced with pitch-maintained slow motion audio, foreign language learning, 
telephone answering machines, and post-production for the film industry. 

TSM algorithms fall into three main categories, time domain approaches, 
frequency domain approaches, and parametric modelling approaches. The 
simplest (and most computationally efficient) algorithms are time domain ones 
and nearly all are based on the principal of Overiap Add (OLA) or Synchronous 
Overiap Add (SOLA), as described in "Non-parametric techniques for pitch scale 
and time scale modification of speech" by E. Moulines and J. Laroche, Speech 
Communications, Vol.16, 1995, pp1 75-205, and "An Edge Detection Method for 
Time Scale Modification of Accoustic Signals" by Rui Ren of the Hong Kong 
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University of Science & Technology Computer Science Department, viewed at 
http://\AAAAA^.cs.ust.hk/-rren/sound tech/TSM_Paper_Long.htm . In OLA, a short 
time frame of music or speech containing several pitch periods of the 
fundamental frequency has a predetermined length: to increase this, a copy of 
5 the input short time frame is overlapped and added to the original, with a cross- 
fade applied across this overlap to remove discontinuities at the block 
boundaries, as will be described in greater detail hereinafter with reference to 
Figures 2, 3 and 4. Although the OLA procedure is simple and efficient to 
implement, the resulting quality is relatively poor because reverberation effects 
10 are introduced at the frame boundaries (splicing points). These artefacts are a 
result of phase information being lost between frames. 
"^-J To overcome these local reverberations, the SOLA technique was 

1=1 proposed by S. Roucos and A. Wilgus in "High Quality Time-Scale Modification 

for Speech", IEEE International Conference on Acoustics, Speech and Signal 
3 15 Processing, March 1985, pp493-496. In this proposal, a rectangular synthesis 

m window was allowed to slide across the analysis window over a restricted range 

2^; generally related to one pitch period of the fundamental. A normalised cross 

p correlation was then used to find the point of maximum similarity between the 

'^^ data blocks. Although the SOLA algorithm produces a perceptually higher 

20 quality output, the computational cost required to implement the normalised 
cross correlation make it impractical for systems where memory and CPU are 
limited. 

It is an object of the present invention to provide a signal processing 
25 technique (and an apparatus employing the same) which, whilst based on 
SOLA techniques, provides a similar quality at a lower computational cost. 

In accordance with the present invention there is provided a method of 
time-scale modification processing of frame-based digital audio signals wherein, 
for each frame of predetermined duration: the original frame of digital audio is 
30 copied; the original and copied frames are partly overiapped to give a desired 
new duration to within a predetermined tolerance; the extent of overiap is 
adjusted within the predetermined tolerance by reference to a cross correlation 
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determination of the best match between the overlapping portions of the original 
and copied frame; and a new audio frame is generated from the non- 
overlapping portions of the original and copied frame and by cross-fading 
between the overlapping portions; 
5 characterised in that a profiling procedure is applied to the overlapping 

portions of the original and copied frame prior to cross correlation, which 
profiling procedure reduces the specification of the respective audio frame 
portions to respective finite arrays of values, and the cross correlation is then 
performed in relation only to the pair of finite arrays of values. By the 

10 introduction of this profiling procedure, the volume of data to be handled by the 
computationally intensive cross correlation is greatly reduced, thereby permitting 
implementation of the technique by systems having lower CPU and/or memory 
capability than has heretofore been the case. 

For the said overlapping portions the profiling procedure suitably identifies 

15 periodic or aperiodic maxima and minima of the audio signal portions and places 
these values in the respective arrays. For further ease of processing, the 
overlapping portions may each be specified in the form of a respective matrix 
having a respective column for each audio sampling period within the 
overlapping portion and a respective row for each discrete signal level specified, 

20 with the cross correlation then being applied to the pair of matrices. A median 
level may be specified for the audio signal level, with said maxima and minima 
being specified as positive or negative values with respect to this median value. 

To reduce computational loading, prior to cross correlation, at least one of 
the matrices may be converted to a one-dimensional vector populated with 

25 zeros except at maxima or minima locations for which it is populated with the 
respective maxima or minima magnitude. 

In the current implementation, the maximum predetermined tolerance 
within which the overlap between the original and copied frames may be 
adjusted suitably, has been restricted to a value based on the pitch period (as 

30 will be described in detail hereinafter) of the audio signal for the original frame to 
avoid excessive delays due to cross correlation. Where the aforesaid median 
value is specified, the maxima or minima may be identified as the greatest 
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recorded magnitude of the signal, positive or negative, between a pair of 
crossing points of said median value: a zero crossing point for said median 
value may be determined to have occurred when there is a change in sign 
between adjacent digital sample values or when a signal sample value exactly 
5 matches said median value. 

Also in accordance with the present invention there is provided a digital 
signal processing apparatus arranged to apply the time scale modification 
processing method recited above to a plurality of frames of stored digital audio 
signals, the apparatus comprising storage means arranged to store said audio 
10 frames and a processor programmed, for each frame, to perform the steps of: 

copying an original frame of digital audio and partly overlapping the 
original and copied frames to give a desired new duration to within a 
predetermined tolerance; 

adjusting the extent of overlap within the predetermined tolerance by 
15 applying a cross correlation to determine the best match between the 
overlapping portions of the original and copied frame; and 

generating a new audio frame from the non-overlapping portions of the 
original and copied frame and by cross-fading between the overlapping portions; 

characterised in that the processor is further programmed to apply a 
20 profiling procedure to the overlapping portions of the original and copied frame 
prior to cross correlation to reduce the specification of the respective audio 
frame portions to respective finite arrays of values, and apply the cross 
correlation in relation only to the pair of finite arrays of values. 

25 Further features and preferred embodiments of the present invention will 

now be described, by way of example only, and with reference to the 

accompanying drawings, in which: 

Figure 1 is a block schematic diagram of a programmable data 

processing apparatus suitable to host the present invention; 
30 Figure 2 illustrates the known Overlap Addition (OLA) time extension 

process; 

Figure 3 illustrates the matching of audio signal segments from a pair of 




5 PHB 34344 US 

overlapping copies of an audio file; 

Figure 4 represents the loss of phase information at the overlap 
boundary for the signal segments of Figure 3; 

Figure 5 represents the generation of a sparse matrix representation of 
5 an audio signal segment for subsequent cross correlation; 

Figure 6 represents overlap addition for a pitch increase; 

Figure 7 illustrates movement of samples for Time Scale Modification 
buffer management; 

Figure 8 is a table of sample values for analysis and synthesis blocks in a 
10 sparse cross correlation; and 

Figure 9 illustrates in tabular form the progress of a further simplified 
cross correlation procedure. 

Figure 1 represents a programmable audio data processing system, such 

15 as a karaoke machine or personal computer. The system comprises a central 
processing unit (CPU) 10 coupled via an address and data bus 12 to random- 
access (RAM) and read-only (ROM) memory devices 14, 16. The capacity of 
these memory devices may be augmented by providing the system with means 
18 to read from additional memory devices, such as a CD-ROM, which reader 

20 18 doubles as a playback deck for audio data storage devices 20. 

Also coupled to the CPU 10 via bus 12 are first and second interface 
stages 22, 24 respectively for data and audio handling. Coupled to the data 
interface 22 are user controls 26 which may range from a few simple controls to 
a keyboard and a cursor control and selection device such as a mouse or 

25 trackball for a PC implementation. Also coupled to the data interface 22 are one 
or more display devices 28 which may range from a simple LED display to a 
display driver and VDU. 

Coupled to the audio interface 24 are first and second audio inputs 30 
which may (as shown) comprise a pair of microphones. Audio output from the 

30 system is via one or more speakers 32 driven by an audio processing stage 
which may be provided as dedicated stage within the audio interface 24 or it 
may be present in the form of a group of functions implemented by the CPU 10; 
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in addition to providing amplification, the audio processing stage is also 
configured to provide a signal processing capability under the control of (or as a 
part of) the CPU 10 to allow the addition of sound treatments such as echo and, 
in particular, extension through TSM processing. 
5 By way of example, it will be useful to initially summarise the basic 

principles of OLA/SOLA with reference to Figures 2, 3 and 4 before moving onto 
a description of the developments and enhancements of the present invention. 

Considering first a short time frame of music or speech containing 
several pitch periods of the fundamental frequency, and let it's length be N 

10 samples. To increase the length from N to A/' (say 1.75A/), a copy of the input 
short time frame (length N) is overlapped and added to the original, starting at a 
point StOL For the example A/' = 1 .75A/, StOI is 0.75A/. This arrangement is 
shown in Figure 2. The shaded region is the overlap between the data blocks 
(length Of) and, as can be seen from the lower trace, a linear cross fade is 

15 applied across this overlap to remove discontinuities at the block boundaries. 

Although the OI_A procedure is simple and efficient to implement, the 
resulting quality is relatively poor because reverberation effects are introduced 
at the frame boundaries (splicing points). These artefacts are a result of phase 
information being lost between frames. 

20 In the region of the overlap we define the following. The analysis block is 

the section of the original frame that is going to be faded out. The synthesis 
block is the section of the overlapping frame that is going to be faded in (i.e. the 
start of the audio frame). The analysis and synthesis blocks are shown in 
Figure 3 at (a) and (b) respectively. As can be seen, both blocks contain similar 

25 pitch information, but the synthesis block is out of phase with the analysis block. 
This leads to reverberation artefacts, as mentioned above, and as shown in 
Figure 4. 

To overcome these local reverberations, the SOLA technique may be 
applied. In this technique, a rectangular synthesis window is allowed to slide 
30 across the analysis window over a restricted range [0, K^aJ where K^^ 
represents one pitch period of the fundamental. A normalised cross correlation 
is then used to find the point of maximum similarity between the data blocks. 
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The result of pitch synchronisation is shown by the clashed plot in Figure 3 at 
(c). The synthesis waveform of (b) has been shifted to the left to align the peaks 
in both waveforms. 

As mentioned previously, although the SOLA algorithm produces a 
perceptually high quality output, the computational cost required to implement 
the normalised cross correlation make it impractical to implement for systems 
where CPU and memory are limited. Accordingly, the present applicants have 
recognised that some means is required for reducing the complexity of the 
process to allow for its implementation in relatively lower powered systems. 

The normalised cross correlation used in the SOLA algorithm has the 
following form: 



Rik) 



^-0,1,2, 



(1) 



2 



V J 



where j is calculated over the range [0, O/], where Ol is the length of the 
15 overlap, x is the analysis block, and y is the synthesis block. The maximum 
R(k) is the synchronisation point. 

In terms of processing, this requires 3x0/ multiply accumulates (macs), 
one multiply, one divide and one square root operation per k value. As the 
maximum overlap that is considered workable is 0.95A/, the procedure can 
20 result in a huge computational load. 

Ideally the range of k should be greater than or equal to one pitch 
period of the lowest frequency that is to be synchronised. The proposed value 
for Kmax in the present case is 448 samples. This gives an equivalent pitch 
synchronising period of approximately 100 Hz. This has been determined 
25 experimentally to result in suitable audio quality for the desired application. For 
this k value, the normalised cross correlation search could require up to 
approximately 3 million macs per frame. The solution to this excessive 
number of operations consists of a profiling stage and a sparse cross 
correlation stage, both of which are discussed below. 
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Both the analysis and synthesis blocks are profiled. This stage consists 
of searching through the data blocks to find zero crossings and returning the 
locations and magnitudes of the local maxima and minima between each pair 
of zero crossings. Each local maxima (or minima) is defined as a profile point. 
5 The search is terminated when either the entire data block has been 
searched, or a maximum number of profile points (P^ax) have been found. 

The profile information for the synthesis vector is then used to generate 
a matrix, S with length equal to the profile block, but with all elements initially 
set to zero. The matrix is then sparsely populated with non-zero entries 
10 corresponding to the profile points. Both the synthesis block 100 and S are 
shown in Figure 5. 

It is clear from this example that the synthesis block has been replaced 
yi by a matrix S which contains only six non-zero entries (profile points) as 

%l shown at 101-106. 

'^i 15 In order to determine the local maxima (or minima) between zero 

p crossings, the conditions for a zero crossing must be clearly defined. 

™J Subjective testing with various configurations of zero crossing have led to the 

nJ following definition of a zero crossing as occurring when there is either: 

p • a change in sign from a positive non-zero number to a negative non-zero 

20 number, and vice versa; or 

• there is an element with a magnitude of exactly zero. 

Transitions from positive to zero or from negative to zero are not included in 

the definition. 

Turning now to calculating the sparse cross correlation, the steps 
25 involved are as follows. Firstly, both the analysis and synthesis waveforms 
are profiled. This results in two 2-D arrays Xp and Yp respectively, of the form 
Xp(loc, mag), where: 

Xp(0,0) = location of first maxima (minima), 
Xp(0,1) = magnitude of first maxima (minima). 
30 Each column of the profiled arrays contains the location of a local 

maxima (or minima) and the magnitude of the maxima (or minima). These 
arrays have length = P^naiysis or Psynthesis^ and a maximum length = P^^,, the 
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maximum number of profile points. 

A 1-D synthesis vector S (which has length = length of synthesis buffer) 
is populated with zeros, except at the locations in yp(/,0), where / = 0,1,... 
"■Psynthesisy where it is populated with the magnitude y(/,1). 

The sparse cross correlation now becomes: 



R\k) = 



Pa na/ys is-l 

^x(/,l)x s{x{i,0)-\- k) 



i=0 



f Panalysis-l \/ Ploc-\ 



(2) 



where P,oc is the number of synthesis points that lie within the range [0+k, 
Ol+k]. 

10 As can be seen, the square root has been removed. Also it can be 

' analysis 

seen that the energy calculation ^jcj only needs to be calculated once a 

J 

frame and so can be removed from equation 2. 

The resulting number of macs required per frame is now limited by the 
maximum number of analysis profile points {PnjaxY a preferred 

15 implementation, P^^^ = 127, which has been found to provide ample resolution 
for the search. This means that for each frame, the Worst Case Computational 
Load per frame = 2x127x448 is limited now by P^rjax » as opposed to OA The 
improvement factor can be approximated by Ol / P^a^ which, for an overlap of 
2048 samples, results in a reduction of the computational load by a factor of 

20 approximately 10. There is an additional load of approximately 12.5k cycles 
per frame, but this is of the order of 20 to 30% improvement in computational 
efficiency. Both objective and informal subjective tests performed on the 
present method and SOLA algorithm produced similar results. 

Considering now the issue of buffer management for the TSM process, 

25 overlapping the frames to within a tolerance of K^ax adds the constraint that 
the synthesis buffer must have length = O/ + K^^x- As this is a real-time 
system, another constraint is that the time scale block must output a minimum 
of A/' samples every frame. To allow for both constraints the following buffer 
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management is implemented. The cases for pitch increases and pitch 
decreases are different and so will be discussed separately. 

Considering pitch increase initially, Figure 6 shows the process of time 
expansion with pitch synchronisation. It is apparent from the diagram that if k 
5 = Krrjax, the length of the time extended frame will be less than N\ To solve 
this, StOI is simply increased by K^^^. This results in spare samples (in the 
range [0, K^axl) the end of the frame. These samples are stored in a buffer 
and added on to the start of the next frame as shown in Figure 7. This results 
in a variable length {N^ctuaD ^or the current input frame, so the scale factor (i.e. 
-ij 10 A/7 Nsctuai) nnust be recalculated every frame. If for a given frame N ever 
Si exceeds N\ then A/' samples from the input frame are outputted and any 

5^ remaining samples are added onto the start of the next frame. 

Q Turning now to pitch decrease, in this case samples remaining from the 

^ ' previous frame are stored and overlapped and added to the start of the current 

^ 15 frame. The analysis block is now the start of the current frame, and the 
nJ synthesis block is comprised of samples from the previous frame. Again, the 

P synthesis block must have length greater than Ol + K^ax'"^- 'f the synthesis 

block is less than this length it is simply added onto the start of the current 
input frame. A/' samples are outputted, and the remaining samples are stored 
20 to be synchronously overlap added to the next frame. This procedure 
guarantees a minimum of A/' samples every frame. 

In order to allow a smooth transition between frames a linear cross fade 
is applied over the overlap. This cross fade has been set with two limits; a 
minimum and a maximum length. The minimum length has been determined 
25 as the length below which the audio quality deteriorates to an unacceptable 
level. The maximum limit has been included to prevent unnecessary load 
being added to the system. In this implementation, the minimum cross fade 
length has been set as 500 samples and the maximum has been set at 1000 
samples. 
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A further simplification that may be applied to improve the efficiency of 
the sparse cross correlation will now be described with reference to the tables 
of Figures 8 and 9. 

Consider first the table of Figure 8 which shows the results of profiling 
5 the analysis and synthesis frames. Arrays Sp and Ap are created (from the 
synthesis and analysis frames respectively), each of which holds a maximum 
of 127 profile entries, each entry containing the magnitude of the profile point, 
as well as the location at which that point was found in the original analysis 
and synthesis frames. This is different from the earlier implementation, in that 
P== 10 only one low entry profile array was created, and the other frame (the 
;fl synthesis frame) was represented by a sparsely populated array of the same 

size as the original frame. As can be seen from the Figure, each array is 

yi 

m terminated with -1 in the location entry to indicate the profile is complete. 

—f 

Si In order to calculate the profile, for each value of j=O..K , the following is 

fi 15 undertaken: 

^5 Initialise variables Ap_count and Sp_count to zero, 

nj Chose either Ap or Sp (say Ap) as the initial driving array. Driving and 

5^ non driving arrays d and nd are provided as pointers, which are then used to 

point to whichever of Ap or Sp are the driver for a particular iteration through 
20 the algorithm. These also hold values d_count and nd_count, which are used 
to hold the intermediate values of ap_count and sp_count whilst a particular 
array is serving as the driving array. 

It will be noted that, depending upon which array is the driving array, in 
practice either the .loc or Joe +k value is used in later calculations. This may 
25 be done efficiently, for example, by always adding j*gate to the .loc value, and 
gate is a value either 0 or 1 depending upon whether the analysis frame is 
chosen. So, d_gate and nd__gate, hold these gate values and when the 
driving array pointer is swapped the gate values should also be swapped 
Hence a comparison of the loc values of the driving and non-driving arrays 
30 will be: 

Is driving[d_count].loc + j*d_gate > non_driving[nd_count].loc + j*nd_gate 
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So, starting to perform an iteration. 

Compare driving[d_count].loc + j*d_gate with non_driving[nd_count].loc + 
j*nd_gate. 

5 

If the two locations match, either perform the cross correlation summations 
now, or else add the Ap and Sp magnitude values (accessed in the same 
manner as the .loc values) to a list of Values to multiply later'. Increment 
Sp_count and Ap_count (d count and nd_count), and pick a new driving array 
10 by finding the maximum of the numbers Ap[Ap__count].loc, Sp[Ap_count].loc + 
j (if the two match then pick either), thus giving a new driving array to guide 
the calculations. 

If the values do not match, then: 

• if the .loc value in the driving array is greater than the .loc value in the non- 
15 driving array, then increment the _count value of the non driving array. 

• If the .loc of the driving array is less than the .loc of the non-driving array 
then increment the _count value of the driving array 

• Make the driving array the one with the higher loc value, unless both are 
the same, in which case do nothing. 

20 Now perform a new iteration and continue with this until either array is —1 
terminated, indicating one of the profile arrays is exhausted. If the 
multiplications were not performed during the above phase, the list of 
magnitude values to multiply together should now be extracted and the cross- 
correlation calculated. In the example above, the process is illustrated for j = 

25 1. 

In the above approach only two multiplications are carried out j=1, as 
compared to a total of 4 which would have been required in a dumb 
implementation, with the added complexity of the implementation above. On 
the face of it this is an insignificant depreciation, but, as the number of profile 
30 points increase, then the scope for reducing the number of multiplications 
decreases further. Effectively the number of multplications that are carried out 
is bounded by the smaller of the number of points in either profile array, as 
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opposed to being bounded by the number in the analysis array as in the 
earlier implementation, which gives potential for high gains. 

Although defined principally in terms of a software implementation, the 
skilled reader will be well aware than many of the above-described functional 
5 features could equally well be implemented in hardware. Although profiling, 
used to speed up the cross correlation, dramatically reduces the number of 
macs required, it introduces a certain amount of pointer arithmetic. Processors 
such as the Philips Semiconductors TriMedia ™ , with its multiple integer and 
floating point execution units, is well suited to implementing this floating point 

10 arithmetic efficiently in conjunction with floating point macs. 

The techniques described herein have further advantage on TriMedia in 
that it makes good use of the TriMedia cache. If a straightforward cross 
correlation were undertaken, with frame sizes of 2*2048, it would require 16k 
data, or a full cache. As a result there is likely to be some unwanted cache 

15 traffic. The approach described herein reduces the amount of data to be 
processed as a first step, thus yielding good cache performance. 

From reading the present disclosure, other modifications will be apparent 
to persons skilled in the art. Such modifications may involve other features 
20 which are already known in the design, manufacture and use of image 
processing and/or data network access apparatus and devices and component 
parts thereof and which may be used instead of or in addition to features 
already described herein. 



